Selection rules in a strongly coupled qubit-resonator system 
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Superconducting qubits acting as artificial two-level atoms allow for controlled variation of the 
symmetry properties which govern the selection rules for single and multiphoton excitation. We 
spectroscopically analyze a superconducting qubit-resonator system in the strong coupling regime 
under one- and two-photon driving. Our results provide clear experimental evidence for the con- 
trolled transition from an operating point governed by dipolar selection rules to a regime where 
one- and two-photon excitations of the artificial atom coexist. We find that the vacuum coupling 
between qubit and resonator can be straightforwardly extracted from the two-photon spectra where 
the detuned two-photon drive does not populate the relevant resonator mode significantly. 



I. INTRODUCTION 



Superconducting quantum circuits have proven to be ver- 
satile model systems for performing quantum optics ex- 
periments on a chip, allowing for the investigation of 
fundamental quantum mechanics and scenarios for scal- 
able quantum information processing. With respect to 
the former, circuit quantum electrodynamics (QED)i*'^ 
has gained exceptional importance in recent years. Here, 
the coupling of natural atoms to the photon field in- 
side a three-dimensional optical cavity, knovi^n as cav- 
ity QEDj-Slis modeled using superconducting circuits. In 
particular, suitably designed Josephson junction based 
quantum circuits act as artificial two-level atoms (qubits) 
and transmission line resonators serve as "boxes" for mi- 
crowave photons. Therefore, circuit QED systems are ex- 
tremely attractive because high qubit-resonator coupling 
strengths can be attained due to huge effective dipole mo- 
ments of the artificial atoms and the small mode volumes 
of the resonators. For this reason, the strong coupling 
regime, where qubit and resonator exchange photons co- 
herently, can be realized easily despite the limited coher- 
ence times typically ranging from a few nano- to a few mi- 
croseconds. Recently, even the new regime of ultrastrong 
light matter coupling has been reached,'-^ where the cou- 
pling strength becomes a significant fraction of the total 
system energy. Another important key feature of circuit 
QED is the design flexibility inherent to circuits on a 
chip and the ability to tune their properties in situ over 
a wide range of parameters. This has triggered exper- 
imental studies of multi-photon driven artificial atoms, 
including population inversion,!^ Mach-Zehnder interfer- 
onretry,^*^ Landau-Zener interference,'^ qubits coupled to 
microscopic defects,!^ amplitude spectroscopy,Q^ multi- 
photon spectroscopy of hybrid quantum systems,131 and 
sideband transitions.^^ Recently, an experimental inves- 
tigation of a two-photon driven flux qubit coupled to a 
lumped-element LC resonator demonstrated the control- 



lability of fu ndam ental symmetry properties of circuit 
QED syst ems .f^SI^ However, in that work, a high loss rate 
of the LC resonator leading to a weakly coupled qubit- 
resonator system in combination with the presence of mi- 
croscopic defects complicated the analysis of the data. 
Here, we report on one- and two-photon spectroscopy 
of a qubit-resonator system in the strong coupling limit. 
Our results provide clear experimental evidence for the 
controlled transition from an operating point governed by 
dipolar selection rule^^^ to a regime where one- and two- 
photon excitations of the artificial atom coexist. While 
the former is strictly analogous to the case of natural 
atoms, the latter regime is a particular prominent fea- 
ture of circuit QED. In addition, we find that the vac- 
uum coupling strength between qubit and resonator can 
be straightforwardly extracted from the two-photon spec- 
tra since in this case the detuned two-photon drive does 
not populate the relevant resonator mode significantly. 
The paper is structured as follows: in Sec. [llj we intro- 
duce the sample design, fabrication details, measurement 
setup and techniques. In Sec. |III[ we present and discuss 
the results of one- and two-photon spectroscopy of the 
strongly coupled qubit-resonator system. In Sec. Ill C 
we discuss the symmetry properties of our artificial atom 
and their implications for multi-photon driving in circuit 
QED. Finally, Sec.llVlis dedicated to concluding remarks. 



II. SETUP AND EXPERIMENTAL 
TECHNIQUES 

Figure IT] shows our quantum circuit and a schematic of 
the measurement setup. In our experiments, a super- 
conducting persistent flux qubit ^^ ^^ is galvanically^ 
connected to the center conductor of a niobium coplanar 
waveguide resonator (Fig.flla-c). The niobium film with 
a thickness of 100 nm is deposited by dc-magnetron 
sputtering on a thermally oxidized (50 nm) silicon 
substrate. The patterning of the Nb films is done by 
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FIG. 1: (Color online) (a) Photograph of the superconduct- 
ing coplanar waveguide (CPW) resonator. The position of the 
qubit is indicated by the red box while the positions of the in- 
and output capacitors are indicated by black boxes, (b) SEM 
micrograph of a coupling capacitor with an estimated capac- 
itance Ck. ~ 10 fF. The two fingers are separated by a gap of 
4 iim and are approximately 100 iim long, (c) False-color SEM 
micrograph of the two flux qubits (red) galvanically connected 
to the center conductor of a niobium CPW resonator (blue) 
fabricated on a silicon substrate (gray). The qubit loop areas 
are approximately 8.5 x 14 /im^. (d) Schematic of the mea- 
surement setup. The coupled qubit-cavity system is probed 
through a highly attenuated input line using a vector net- 
work analyzer (VNA). The transmitted signal is amplified at 
4K with a cryogenic HEMT amplifier. A second microwave 
source excites the system in spectroscopy experiments. 



a current antinode of the A-mode (a;2/27r = 5.324 GHz), 
a 80 iim section of the center conductor is replaced by a 
narrow (500 nm) aluminum strip. Two flux qubits are 
galvanically coupled to this narrow strip (see Fig. fl]c). 
The inhomogeneous geometry of the transmission line 
leads to a position dependent inductance L{x) and 
capacitance C{x) per unit length which in turn results 
in a non-harmonic distribution of the higher mode 
frequencies, ujj ^ j ■ loi with j £ N. Each flux qubit 
consists of three Al/AlOx/Al Josephson junctions (JJ) 
interrupting a superconducting loop. The qubits and 
the narrow aluminum strip are fabricated by electron 
beam lithography and shadow evaporation technique. 
The bottom and top aluminum layer have a thickness 
of 50 nm and 70 nm, respectively. Two of the J J have 
nominal areas A sa 0.03 /im^, while one junction is 
smaller by a factor of a « 0.63. A schematic of the 
measurement setup is shown in Fig. [l]d. The sample is 
placed at the base temperature of 15 mK in a dilution re- 
frigerator. We monitor the amplified cavity transmission 
signal at Wj-f of the coupled qubit-resonator system using 
a vector network analyzer (VNA). For spectroscopy 
measurements, the qubit can be excited by an additional 
microwave tone cOs with power Pg. A cryogenic circulator 
prevents amplifier noise from entering our experiment. 

The potential of the three-junction flux qubit can be re- 
duced to a one-dimensional double well potential, where 
the two minima correspond to opposite persistent cur- 
rents ±/p. This double well potential can be tuned by 
an external flux bias <i>x. At (5<i>x = <i>x - (^ - 0.5)$o = 0, 
these current states are degenerate. Here, $o = h/2e 
denotes the magnetic flux quantum and fc e Z. The 
finite tunnel coupling between these states results in a 
symmetric and antisymmetric superposition state sep- 
arated by an energy gap A. The two flux qubits are 
fabricated with slightly different loop areas by varying 
the electron exposure dose during electron beam lithog- 
raphy. For |fc| > 2, the small area difference allows to 
investigate the individual qubit parameters separately by 
making use of the <i>o-periodicity of the qubit potential. 
The spectroscopy measurements presented in this work 
are conducted around <I>x ~ —2.5 $o {k = —2), where 
the qubit transmission spectra are well separated in the 
relevant experimental frequency regime. In the remain- 
der of this paper we solely focus on one of the qubits,'211 
which we refer to as qubit 1. The qubit-resonator system 
can be described quantum mechanically by a multi-mode 
Jaynes-Cummings Hamiltonian 



optical lithography and reactive ion etching. The center 
conductor has a width of 20 pm and is separated from 
the lateral ground planes by a 12 jim gap, leading to 
a characteristic impedance of about 50 il. The high-Q 
resonant cavity (3500 ^ Q < 6300; depending on the 
mode frequency) is coupled via small gap capacitors C^ 
to the readout circuitry and has a fundamental reso- 
nance frequency (A/2-mode) oi Ldi/2Tr = 2.745 GHz. At 
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Here, Wq = y^A^ -|- (2/p • (5(f>x)^/^ is the flux dependent 
qubit transition frequency and a^ a Pauli spin operator. 
The operator cv, (gj) creates (annihilates) a photon in 
the resonator mode with frequency ujj while a^ (f''-) 
creates (annihilates) an excitation in the qubit. The 
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FIG. 2: (Color online) (a) Sketch of the current distribution 
of the first three modes of the coplanar waveguide resonator 
(CPW) with length L = 23 mm. The resonance frequencies 
at $x = are: uii/2-k = 2.745 GHz (A/2-mode, dashed red), 
aj2/27r = 5.324 GHz (A-mode, blue) and aj3/27r = 7.775 GHz 
(3A/2-mode, green). The current distribution for the A/2- 
mode is plotted dashed as its presence can be safely neglected 
for the qubit under consideration, (b) Cavity transmission 
spectrum (linear scale, arb. units) for the 3A/2-mode (top 
panel) and the A-mode (bottom panel) as a function of the 
relative flux bias (5"I>x of qubit 1 and the probe frequency uj^i . 
The spectra are recorded at a flux bias of $x ~ —2.5 $0, 
where the two qubits present in the system can be studied 
independently. The solid lines show the numerically evalu- 
ated full energy level spectrum, where we included the second 
qubit into Hamiltonian ([l]) to verify our assumption of neg- 
ligible qubit-qubit coupling.'^ In contrast to the black lines, 
the gray lines represent energy levels of states with more than 
one excitation .12 



green line in Fig. [2]a, respectively). For the qubit under 
consideration, it is perfectly fine to neglect the presence 
of the A/2-niode (j = 1; dashed red line in Fig. [2]a) in 
the theoretical modeling as it is both largely detuned 
from Wq and only weakly coupled. The same is true for 
the J > 4 modes. Figure [2]b shows the cavity trans- 
mission spectra for the A-and 3A/2-mode as a function 
of (5$x. The transmission spectrum of the 3A/2-mode 
(Fig. [2]b, top panel) reveals avoided crossings at flux po- 
sitions (<5<i>x ~ i2 m$o) where the transition frequency of 
the qubit Wq equals the resonator mode frequency uj^. At 
these flux values, the system is resonant and energy can 
be coherently exchanged between qubit and resonator. 
There are no avoided crossings in the transmission spec- 
trum of the A-mode for qubit 1 (Fig.[2]b, bottom panel). 
Here, the system is in the dispersive regime for all (5$x 
values. In this regime the detuning 5j = jwq — ujj\ be- 
tween the qubit and the jth cavity mode is large in 
comparison to gj. Although qubit and resonator can 
not exchange energy directly, they are still interacting 
through a second order dispersive coupling. Therefore, 
we observe an ac-Zeeman shift of the cavity's resonance 
frequency W2 which depends on the qubit state. Thus 
the ac-Zeeman shift ca n be utilized for two-tone spec- 
troscopy measurements.'SES The measurement protocol 
for this spectroscopy technique is as follows: for flxed 
(5$x, we monitor both phase and amplitude at the probe 
tone frequency Wrf where maximum cavity transmission 
occurs. A second microwave tone, the spectroscopy tone 
(wg), is applied to the system. When Wg ~ Wq the 
qubit is excited and for sufficiently large amplitude of 
the continuous spectroscopy tone the transition is satu- 
rated, yielding an equal population of ground and excited 
state. This leads to a shift of the cavity's resonance fre- 
quency by gi /6j which results in a decrease in magnitude 
and a phase shift of the transmitted signal at Wrf ■ Fur- 
thermore, the dispersive qubit-cavity interaction modifies 
the qubit 's transition frequency which is then given by 
Wq = Wq + 2g'jajjdj/6j + g'j/Sj^ The last term repre- 
sents the Lamb shift due to the presence of the vacuum. 
The ac-Zeeman shift proportional to fij — {aMj) allows 
to study the influence of cavity photons on the qubit. 
Furthermore, for known gj, a measurement of Wq as a 
function of the probe tone power P^f yields a calibration 
of the intracavity photon number fij using the relation 
Pj.f = fijhujjKj, where Kj/2'k is the full- width at half- 
maximum of the cavity resonance. 



qubit-resonator coupling strength gj = gj sin 9 is flux 
dependent and exhibits a maximum at 6^^ = where 
sin 6 — A/hujq — 1. The maximum qubit-resonator cou- 
pling strength can be expressed as hgj = Mlplj, where 
M denotes the mutual inductance between qubit and cav- 
ity and Ij is the vacuum current of the jth resonator 
mode. Owing to our sample design, gj is maximum for 
the A-mode (see Fig. l2]a). In the sum of the Hamilto- 
nian (n]) we only have to include the experimentally rele- 
vant cavity modes (j = 2, 3; A- and 3A/2-mode, blue and 



III. SPECTROSCOPY IN THE STRONG 
COUPLING LIMIT 

Throughout this section we investigate our coupled 
qubit-cavity system using the A-mode as a readout mode 
for microwave spectroscopy. Figure [Sla shows the dressed 
state level diagram for the single-mode Jaynes-Cummings 
Hamiltonian combining the qubit and the 3A/2-mode 
(a;3/27r = 7.775 GHz). In Fig. [3]b and c, the disper- 
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FIG. 3: (Color online) (a) Dressed state level diagram of the 
coupled qubit-cavity system. The states |±,n) denote the 
eigenstates of a single-mode Jaynes-Cummings Hamiltonian 
with n excitations. The transitions driven in the spectroscopy 
experiments neglecting selection rules are indicated by colored 
arrows: one-photon (magenta), two-photon (green) and two- 
photon blue sideband (bs, blue). A detailed discussion of the 
selec tion r ules applicable to these transitions is given in sec- 
(b) Dressed state level diagram in the dispersive 
> <?3), where |-,1) ^ \g,l) and \ + ,l) -^ |e,0). 
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(c) Dressed state level diagram for the resonant case (^3 — 0). 
Here, the eigenstates |±, 1) are symmetric (-I-) and antisym- 
metric (-) superpositions of \g, 1) and |e,0). 



sive and resonant limits are illustrated in more detail, 
respectively. The arrows in Fig. [3] indicate the one- and 
two-photon transitions which we are driving in the 

spectroscopy experiments neglecting any selection rules. 
For the analysis presented in this section, we introduce 
the notation |±, n) for the eigenstates of the single-mode 
Jaynes-Cummings Hamiltonian, where 



-,n, 



= cosQ\g, n) — sinO|e,ri — 1) 



|-t-, n) = sin0|(7, n) + cos0|e, n — 1) 



(2) 



(3) 



Here, |g,fT.) denotes a state with n photons in the 3A/2- 
mode and the q = {g,e\ describes the qubit ground or 
excited states. The mixing angle Q is given by 



9 = - arctan(2(7j-yn/<5j) 



A. One-photon spectroscopy and determination of 
coupling strengths 

The one-photon spectroscopy data is shown in Fig. |4] and 
reflects the hyperbolic flux dependence of the qubit tran- 
sition frequency. From the spectrum we extract the qubit 
parameters A/ft, = 6.88 GHz and 2/p = 532 nA. Know- 
ing (Dq, the qubit-resonator coupling rates Qj can be de- 
termined by fitting the cavity transmission spectra (see 
Fig. [2]b to the energy level spectrum of the Hamiltonian 
0. The fitted coupfing rates are g2/2TT = 106.5 MHz 
and g-i/'i'K = 90.7 MHz. In the low power limit (Prf, 
Ps -^ 0), the full- width at half maximum (FWHM) of 
the qubit spectroscopy signal is 7 w 50 MHz and all decay 
rates Kj of the resonator modes are smaller than 1.4 MHz. 
Thus, since gj > Kj,j the qubit-resonator system is in 
the strong coupling limit. In order to ensure a correct de- 
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FIG. 4: (Color online) One-photon qubit spectroscopy show- 
ing the hyperbolic flux dependence of Cjq. In the top panel, 
the relative change in the transmission magnitude Sm (color- 
coded) is plotted as a function of (5<l>x and the spectroscopy 
frequency losJI-k. A blue sideband transition (resonance con- 
dition: 2a;s = ixJq + ij;3) is visible. The magenta and blue 
circles in the bottom panel indicate the fitted center frequen- 
cies of the qubit spectroscopy signal and the blue sideband 
spectroscopy signal, respectively. The solid magenta line is 
a fit to the energy level spectrum of Hamiltonian ([l| while 
the solid blue line is evaluated using the fit parameters de- 
rived from the qubit data. The flux-independent feature in 
the region lo^I^-k ~ 7 GHz is attributed to the presence of a 
spurious fluctuator. The data presented in Fig. [2] and Fig. H] 
is taken at an input power of Prf = —134 dBm corresponding 
to n-2. < 0.94. 



termination of A and gj , the transmission spectra and the 
spectroscopy data shown in Fig. [2]b and Fig. |4J respec- 
tively, were recorded in the low-po wer hmit ifij < 0.94). 
The photon number was calibratecP'221 ^y nieasuring the 
photon number dependent ac-Zeeman shift of Wq as de- 
scribed in section HIl Two additional features can be 
clearly identified in the one-photon spectroscopy data: 
first, around ujs/2n sa 7.00 GHz, a flux-independent sig- 
nature is visible and can be attributed to a spurious fluc- 
tuator present in our system. Such fluctuators can be 
either resonant modes or microscopic two-level systems 
and can significantly degrade the coherence properties of 
solid-state based artificial atoms.'22l Furthermore, it was 
shown that the presence of two-level fluctuators can lead 
to symmetry breaking,^^ however, this was not the case 
in our experiment and for the remainder of this paper 
we will therefore neglect the presence of the fiuctuator. 
Second, in the region from 7.4 GHz < uJs/2tt < 7.6 GHz 
our data reveals a flux-dependent spectroscopic signa- 
ture which is consistent with a two-photon driven blue- 
sideband transition of the qubit and the 3A/2-mode. The 
selection rules for such transitions will be discussed in 
section IHICl 



B. Two-photon spectroscopy 

We now investigate the spectroscopic response of the 
qubit-cavity system under direct two-photon driving. 
Again we use the flux dependent transmission maxi- 
mum of the A-mode as probe frequency. However, for 
a two-photon drive, the spectroscopy tone is applied in 
a frequency range around 2ws = Wq. The color-coded 
data and the fitted center frequencies of the spectroscopy 
transmission signal are shown in Fig. [5]a and b, respec- 
tively. The spectrum shows a signature with a hyperbolic 
flux dependence, two anticrossings located symmetrically 
around (5$x = 0, and additional flux-independent fea- 
tures in close vicinity of the anticrossing regions. A flt of 
the spectrum (see Fig. Islb) yields the two-photon qubit 
parameters A^p^/Zi = 3.48 GHz and 2/2ph = 296 nA 
which are in good agreement with the one-photon spec- 
troscopy data. We attribute the deviations from the ex- 
pected values (A^Pi^ = A/2 and /^p^^ = Ip/2) to a higher 
intracavity photon number (^2 ~ 2.2) compared to the 
one-photon spectroscopy measurement, thereby inducing 
an ac-Zeeman shift in Cjq towards higher frequencies. 



C. Symmetry properties and selection rules 

The gradual disappearance of the qubit spectroscopy sig- 
nal (see Fig. [sja) for S^^ — >■ can be understood by con- 
sidering the symmetry properties of our system which im- 
ply selection rules for the allowed transitions. Selection 
rules are intimately related to the quantum mechanical 
concept of parity. The parity operator H divides the set 
of all states into three groups: even states (|V'-i-)), odd 



states (IV'-)) a-nd states without well-defined parity. The 
former two are eigenstates of H with eigenvalues -1-1 and 
— 1, respectively. In a similar fashion, it is possible to 
divide the set of all operators into similar classes. While 
an even operator A^ commutes with H, an odd opera- 
tor A_ anticommutes with H. It can be shown,ES that 
the matrix elements of an even operator are zero between 
states of different parity 



(^+|i+|V'_) = (^_|i+|^+)=0 



(5) 



while the matrix elements of an odd operator are zero 
between states of the same parity 



(V+|^-|^+) = (^-|A_|v_) = o. 



(6) 



At S^x — 0, the two lowest qubit energy eigenstates \g) 
and |e) represent symmetric and antisymmetric super- 
positions of the persistent current states | ± /p). Tran- 
sitions between \g) and |e) are induced by an external 
drive oc (a' -I- a) with respect to a frame rotating at the 
drive frequency. It is easy to verify that the drive {a^ +d) 
is an odd parity operator and as \g) and |e) correspond 
to states with different parities, one-photon transitions 
are allowed at the qubit degeneracy point. On the other 
hand, a two-photon (n-photon) transition is equivalent 
to the application of two (n) subsequent drives with fre- 
quency (jJq/2 {ujq/n). With 



(fli_n)" = (-l)"i!! 



(7) 



the two-photon drive is effectively an even parity op- 
erator and two-photon transitions between \g) and |e) 
are forbidden. Away from the qubit degeneracy point 
((5<I>x 7^ 0), the parity of the Hamiltonian (II]) is not well- 
defined. In the dispersive limit, the effective Hamiltonian 
can be derived^^ "^^ by means of a Schrieffer- Wolff trans- 
formation yielding finite transition matrix elements for 
one photon 



Vt 



and two photon drive 



n^ 



4A 



■ smf 



sin 6 cos 6 



(8) 



(9) 



When deriving the above expressions, we assume the 
drive Hamiltonian H^ = {il/2)azC0s{ujt). Here, w = 
LOq/n for an n-photon drive and Q refers to the driv- 
ing strength. The cos ^-dependence of the two-photon 
spectroscopy signal as qualitatively reproduced by our 
data (see Fig. Isla) explains the gradual disappearance 
of the qubit signature as cos 6 — >• for S^^ — ?> 0. The 
selection rules for one- and two-photon driven red and 
blue sideband transitions can be derived in a similar 
fashion. ^^ The two-photon blue sideband illustrated in 
Fig. I4] corresponds to the transition \g,0) — ^ |e, 1). These 
states have the same parity and since the two-photon 
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FIG. 5; (Color online) (a) Two-photon spectroscopy. The relative change in the transmission magnitude at the probe tone 
frequency ujif is plotted as a function of 5<l>x and the spectroscopy frequency ujs/2n. At 5$x = 0, two-photon transitions are 
strictly forbidden due to selection rules originating in the symmetry of the qubit potential. Away from the degeneracy point, the 
magnitude of the spectroscopy signal increases. The two anticrossings are the spectroscopic correspondents of the qubit-cavity 
anticrossings visible in the cavity transmission spectrum of the 3A/2-mode (see Fig. [2]b, top panel). On the left hand side, a 
one-photon driven red sideband transition (resonance condition: Ws — ujq ~ UJ2) is visible. This transition corresponds to an 
exchange of an excitation between the qubit and the A-mode. (b) Fitted center frequencies of the spectroscopic qubit response 
(green circles) as a function of 5$x. The solid green line is the numerically evaluated energy level spectrum of Hamiltonian (IT]) 
with Ldq — >■ u}ci/2,u)n — > ^n/2 and gn — ^ g-nfi. Dotted blue line: 1/2 • uj3/2tt. Dashed red line: calculated flux dependence of 
the red sideband transition, (c) Two-photon spectroscopy signal in the anticrossing region plotted as a function of u}b/2tt. The 
curves for different constant (S$x are shifted along the ordinate. The gray dashed box marks the area magnified in (d). (d) 
Two-photon spectroscopy signal at a flux value where the qubit is resonant with the 3A/2-mode. The red solid line represents 
a fit to the sum of two Lorentzians from which we extract a dip separation of §3 = 77.4 MHz. 



drive has even parity, the transition is allowed for ar- 
bitrary (S$x- However, because this is a second-order 
transition, its amplitude is small and the spectroscopic 
signature can not be resolved within our measurement 
resolution and for low spectroscopy power Pg. Finally 
we note that for the harmonic resonator potential multi- 
photon transitions are forbidden, ^^ which is illustrated by 
the absence of a flux-independent spectroscopic feature 
at Wg = W3/2 except for the anticrossing regions. Here, a 



flux-independent resonator-like signature is visible due to 
the large qubit-like component of the eigenstatcs of the 
Jaynes-Cummings Hamiltonian (see Eqs. (pi) and (Is])). 



D. Anticrossing under two-photon driving 

The anticrossings in the spectrum shown in Fig. [5]a are a 
direct manifestation of the one-photon qubit-resonator 



anticrossing (see Fig. [2]a) under two-photon driving. 
Here, the qubit and resonator energies are degenerate, 
and the eigenstates |±, 1) of the coupled system are sym- 
metric and antisymmetric superpositions of \g, 1) and 
|e, 0) (see Fig.lslc). The fact that these anticrossings are 
resolved so clearly indicates that the qubit-cavity system 
is in the strong coupling limit. There are several ways 
to determine gj, all requiring a negligible photon pop- 
ulation of the cavity modes. A straightforward way to 
investigate the coupling is to record the vacuum Rabi 
splitting. Observing this splitting is experimentally chal- 
lenging but has been demonstrated with different kinds of 
qubits. '^''^'^^"^^' Typically, the amplitude of the two peaks 
corresponding to the entangled qubit-cavity states is sig- 
nificantly reduced compared to the bare cavity mode due 
to the following reason: on resonance, the two systems 
can exchange energy and a decay of the qubit into non- 
radiative channels eliminates the shared excitation of the 
joint system. Such processes thus lead to a significant re- 
duction of the transmitted microwave amplitude which 
in turn increases the required signal averaging times. 
Furthermore, the Rabi peaks become inhomogeneously 
broadened with increasing power. The coupling can also 
be extracted from a fit of the energy level spectrum of 
the qubit-cavity Hamiltonian (see section III A) to the 



measured cavity transmission spectra. Here, we extract 
the qubit-cavity coupling from the dispersive two-photon 
spectroscopy data. This has the advantage that both 
drive frequencies, uj-^i and uj^ are applied off-resonant to 
the cavity mode resonant with the qubit (in our case 
the 3A/2-mode). Therefore, this mode is not populated 
significantly by the drives and the vacuum Rabi splitting 
can be determined form the spectroscopy data. Figure[5]c 
shows the measured relative change of transmission mag- 
nitude as a function of (5$x and Wg for the anticrossing 
located at (5$x ~ 1.9 m$o- For increasing flux detuning, 
Wq is shifted towards higher frequencies until the reso- 
nance condition Wq = LOy, is fulfilled. We extract a mini- 
mum separation of the two dips of g^ll-n ~ 77.4 MHz as 
shown in Fig. [5]d. Taking into account that sin0 w 0.89 



at this flux value, we obtain (73/27r w 87.0 MHz in good 
agreement with the fitted coupling rate determined from 
the cavity transmission and the one-photon qubit spec- 
troscopy experiments. From the FWHM « (7 -|- k^)/2 of 
the spectroscopy signal we can estimate 7/27r « 43 MHz 
in agreement with the linewidth obtained from low-power 
one-photon spectroscopy. 



IV. CONCLUSIONS 

In summary, we investigate experimentally a super- 
conducting flux qubit strongly coupled to a microwave 
transmission line resonator under one- and two-photon 
driving. At the qubit degeneracy point, the selection 
rules for dipole transitions are equivalent to the selection 
rules for transitions in natural atoms. By tuning 
the external control parameter 5$x the symmetry of 
the system can be broken in a controlled way and 
two-photon driving of the coupled qubit-cavity system 
becomes possible. This technique allows to extract the 
coupling rates g^ of modes resonant with the qubit in 
the dispersive limit. 
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